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Abstract 

In a recent work an approximation procedure was introduced to calculate the vac- 
uum expectation value of the stress-energy tensor for a conformal massless scalar field 
in the classical background determined by a particular colliding plane wave space-time. 
This approximation procedure consists in appropriately modifying the space-time ge- 
ometry throughout the causal past of the collision center. This modification in the 
geometry allows to simplify the boundary conditions involved in the calculation of the 
Hadamard function for the quantum state which represents the vacuum in the flat 
region before the arrival of the waves. In the present work this approximation proce- 
dure is applied to the non-singular Bel-Szekeres solution, which describes the head on 
collision of two electromagnetic plane waves. It is shown that the stress-energy ten- 
sor is unbounded as the killing-Cauchy horizon of the interaction is approached and 
its behavior coincides with a previous calculation in another example of non-singular 
colliding plane wave space-time. 
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1 Introduction 



It is known that exact gravitational plane waves are very simple time dependent plane 
symmetric solutions of Einstein's equations |T|. Nevertheless, they show two main nontrivial 
global features, namely: i) the absence of a global Cauchy surface, which is a consequence of 
the focusing effect that the waves exert on null rays ||, ii) the presence of a Killing-Cauchy 
horizon which may be physically understood as the caustic produced by the focusing of null 
rays [0]. The inverse of the focusing time is a measure of the strength of the wave. For 
an Einstein-Maxwell plane wave such inverse time equals the electromagnetic energy per 
unit surface of the wave. This makes exact plane waves very different from their linearized 
counterparts, which have no focusing points and admit a globally hyperbolic space-time 
structure. One expects that exact plane waves may be relevant for the study of the strong 
time dependent gravitational fields that may be produced in the collision of black holes J3], [5| 
or to represent travelling waves on strongly gravitating cosmic strings ||. In recent years 
these waves have been used in classical general relativity to test some conjectures on the 
stability of Cauchy horizons R M, and in string theory to test classical and quantum string 



behaviour in strong gravitational fields P, 10, O. Their interest also stems from the fact 



that plane waves are a subclass of exact classical solutions to string theory [12, 13, 14 



In Einstein-Maxwell theory the particular class of plane symmetric waves are seen to 
contain only a non-null component of the Ricci tensor and only a non-null component of the 
Weyl tensor. In particular, the single component of the Weyl tensor may be conveniently 
interpreted as the transverse wave component in the direction of propagation of the wave. In 
that sense, the modulus term of the Weyl component can be identified with the amplitude of 
the wave and the phase term with the polarization of the wave. Furthermore, depending on 
whether the Ricci component or the Weyl component is zero we will distinguish in between 
pure gravitational plane waves or pure electromagnetic plane waves respectively. 

When we consider a plane wave collision, we should analyze separately the collision 
between pure gravitational waves, between pure electromagnetic waves or between mixed 
waves. Namely: i) when two pure gravitational plane waves interact, the focusing effect of 
each wave distorts the causal structure of the space-time near the null horizons that these 
waves contain and either a spacelike curvature singularity or a new regular Killing-Cauchy 
horizon is created, ii) when two pure electromagnetic plane waves interact, the situation 
is more subtle. In fact, in the full Einstein- Maxwell theory, Maxwell's equations remain 
linear indicating non direct electromagnetic interaction between the waves. However, there 
is a non-linear interaction of the waves through the gravitational field generated by their 
electromagnetic energy, which is similar to the magnitude of the interaction between pure 
gravitational waves. In that sense, the collision of two electromagenic waves is seen to 
produce gravitational waves, iii) in the case of mixed collisions, the pure electromagnetic 
wave is partially reflected by the incident pure gravitational wave. The gravitational wave, 
however, is not necessarily reflected. 

Note that, the presence of a Killing-Cauchy horizon in a colliding plane wave space-time 
implies a breakdown of predictability since the geometry beyond the horizon is not uniquely 
determined by the initial data posed by the incoming colliding waves. Also, the singularities 
derived from plane wave collisions are not the result of the collapse of matter but the result 
of the non-linear effects of pure gravity, i.e. the self- gravitation of the gravitation field. 

When waves are coupled to quantum fields there is neither vacuum polarization nor the 
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spontaneous creation of particles. In that sense they behave very much as electromagnetic 
or Yang-Mills plane waves in flat space-time |T3|, IB]. However, as a result of the non-linear 
interaction, the creation of quantum particles is expected in a plane wave collision. 

The interaction of quantum fields with colliding plane waves was first considered by 
Yurtsever |T7| for the singular Khan-Penrose solution |I5| , which describes the collision 
of two plane impulsive gravitational waves. In that case, an unambiguous "out" vacuum 
state was possible to define in a relatively simple way. More recently, Dorca and Verdaguer 
H, |(J noticed that the presence of a Killing- Cauchy horizon in a non singular colliding 
plane wave space-time could be used to define an unambiguous "out" vacuum state related 
to the preferred Hadamard state introduced by Kay and Wald in more generic space-times 
with Killing- Cauchy horizons pll . With this premise, Dorca and Verdaguer studied the 
interaction of quantum fields in a particular non-singular colliding plane wave space-time, 
the interaction region of which was isometric to a region inside the event horizon of a 
Schwarzschild black hole p2L |23|. Later on, the same premise was applied by Feinstein and 



Sebastian [24] to the Bel-Szekeres solution |25], which represents the head on collision of two 
electromagnetic plane waves with an interaction region isometric to the Bertotti-Robinson 
universe [26] filled with an uniform electric field. In all these examples it was found that 
the initial state, defined to be the vacuum state in the flat region before the arrival of the 
waves, contained a spectrum of "out" particles consistent, in the long wavelength limit, with 
a thermal spectrum with a temperature inversely proportional to the focusing time of the 
waves. 

A further step in the study of the interaction of quantum fields with colliding plane waves 
is the computation of the expectation value of the stress-energy tensor. Again, this problem 
was first considered by Yurtsever |17| for the Khan-Penrose solution [|Tj|. In that case it 



was possible to determine the behavior of the stress-energy tensor near the singularity of 
the interaction region. It was shown that for the conformal coupling case (i.e. £ = 1/6) the 
energy density and two of the principal pressures were positive and unbounded towards the 
singularity. This problem has been also considered by Dorca and Verdaguer in the mentioned 
above non-singular colliding plane wave spacetime with an interaction region isometric to 
an interior region of a Schwarzschild black hole. As in the case of Yurtsever for the Khan- 
Penrose solution, the expectation value of the stress-energy tensor was calculated in the 
state representing the Minkowski vacuum in the flat region before the arrival of the waves. 
This value was first computed in a region close to both the Killing-Cauchy horizon and 
the topological singularities, the folding singularities, that the colliding plane wave space- 
time contains [28 . In that particular region, the calculations were simplified due to the 



blueshift effect on the energy of the initial quantum modes as they reached the Killing- 
Cauchy horizon ||29|| . In a recent work [[3(J, an approximation procedure was proposed by 
the author in order to calculate such an expectation value throughout the causal past region 
of the collision center. In both calculations, it was found that the stress-energy diverged 
as the Killing-Cauchy horizon was approached. The rest energy density was positive and 
unbounded towards the horizon. Two of the principal pressures were negative and of the 
same order of magnitude of the energy density. It was also pointed out that such a behavior 
suggested that the non singular Killing-Cauchy horizon is indeed unstable under quantum 
perturbations and a curvature singularity would be the general outcome of a generic plane 
wave space-time when backreaction is taking into account. 
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In the present paper, the approximation procedure introduced in [[Kj is applied to the 
non singular Bel-Szekeres space-time as a first attempt to generalize such an approximation 
to more generic colliding plane wave space-times. The plan of the paper is the following. 
In section 2 the geometry of the Bel-Szekeres solution is briefly reviewed. In section 3 an 
adequate approximation in the space-time geometry is introduced throughout the causal past 
of the collision center. Then, the mode solutions of a massless scalar field which represent 
the vacuum state before the arrival of the waves are calculated all over this particular region. 
In section 4 the Hadamard function, which is the key ingredient for the computation of the 
stress-energy tensor, is calculated and regularized by means of the point-splitting technique. 
In section 5 the vacuum expectation value of the stress-energy tensor is calculated. In section 
6 a summary and some consequences of that result are given. In order to help to maintain 
the main body of the paper reasonably clear, some final results are stored in the Appendices. 



2 Description of the geometry 



The Bel-Szekeres solution [25] represents the collision of two electromagnetic shock waves 
followed by trailing radiation. The interaction region is isometric to the Bertotti-Robinson 
universe |p6|| , which is the static conformaly flat solution of Einstein-Maxwell equations with 
an uniform electric field. Such a geometry is similar to the throat of the Reissner-Nordstrom 
solution for the special case M = Q The space-time contains four space-time regions, 
given by 

ds 2 IV = ALiL 2 dudv — dx 2 — dy 2 , (1) 
ds 2 jj = AL\L 2 dudv — cos 2 v [dx 2 + dy 2 ^j , (2) 

dsjj = ALiL 2 dudv — cos 2 u (dx 2 + dy 2 ^j , (3) 
ds 2 = 4LiL 2 dudv — cos 2 (m + v ) dx 2 — cos 2 (m — v) dy 2 , (4) 

where for convenience we have used u and v as dimensionless null coordinates, and where L\ 
and L 2 , are length parameters such that LiL 2 is directly related to the focusing time of the 
collision, i.e. to the inverse of the strength of the waves, which is a measure of the amount 



of nonlinearity of the gravitational waves [p8 



This colliding wave space-time, as shown in Fig. 1, consists of two approaching waves, 
regions II and III, in a flat background, region IV, and an interaction region, region I. 
The two waves move in the direction of two null coordinates u and v, and since they have 
translational symmetry along the transversal x-y planes, the interaction region retains a 
two-parameter symmetry group of motions generated by the Killing vectors d x and d y . The 
four space-time regions are separated by the two null wave fronts u = and v = 0. Namely, 
the boundary between regions I and II is {0 < u < tt/2, v = 0}, the boundary between 
regions I and III is {u — 0, < v < vr/2}, and the boundary of regions II and III with 
region IV is {u < 0, v = 0} U {u = 0, v < 0}. Region I meets region IV only at the surface 
u — v — 0. The Killing-Cauchy horizon in the region I corresponds to the hypersurface 
u + v = 7r/2 and plane wave regions II and III meet such a Killing-Cauchy horizon only at 
V = {u = 7r/2, v = 0} and V = {u = 0, v = vr/2} respectively. Observe that plane wave 
regions II and III contain a singularity at u = tt/2, for region II, and v = tt/2, for region 



3 



III. These singularities are not curvature singularities but a type of topological singularity 
commonly referred to as a folding singularity This terminology arises from the fact 

that the whole singularity u = n/2 in region II (or v = tt/2 in region III) must be identified 
"folded") with V (or V) (see |1| for more details and for a 3-dimensional plot of a 



i.e. 



space-time of this type). 



3 Mode propagation 

For simplicity we will consider in this section a massless scalar field, which satisfies the usual 
Klein-Gordon equation, 

□0 = 0. (5) 
Following the directions of the approximation procedure introduced in the previous work 



30[ , we will be interested in the value of the quantum field all over the causal past 
region of the collision center. The reason is essentially because the calculations can be 
greatly simplified in this region. We will start with the field solution in the flat region prior 
to the arrival of the waves, which is chosen to be the usual vacuum in Minkowski space- 
time. This vacuum solution will set a well posed initial value problem on the null boundary 
£ = {u = 0, v < 0} U {u < 0, v = 0}, by means of which a unique solution for the field 
equation can be found throughout the space-time, i.e., in the plane wave regions (regions II 
and III), and in the interaction region (region I). 
We will consider the line element, 

ds 2 = 2e~ M ^dudv - e~ u ^ (e v ^dx 2 + e- v ^dy 2 ) , (6) 

which applies globally to the four space-time regions, and where the functions U, V and 
M, can be directly read off ([TJ)-(|]). Then, the field equation can be separated in a plane- 
wave form solution for the transversal coordinates x and y, with k x and k y , respectively, 
as separation constants. This plane- wave separation is just a trivial consequence of the 
translational symmetry of the space-time on the planes spanned by the Killing vectors d x 
and d y . The field solution is thus, 

(f>{u, v, x, y) = e u{u ' v)/2 f(u, v) e ik * x+ik * y , (7) 
where the function f(u,v) satisfies the following second order differential equation, 

(e- u / 2 ) i 

f, uv + v) f = 0; n(u, v) = - V e _ uf J 2 uv + l -e- M+u (k 2 x e~ v + k 2 y e v ) . (8) 

Equation (R) can be straightforwardly solved in the flat region (region IV). Then, this 
solution determines on the null boundary S a well posed set of initial conditions for the 
solutions of equation (|§|) in plane wave regions II and III. Finally, the field solution in 
regions II, III and IV is, 
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(u,v,x,y) 
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i2k+ ta.nu—i2k—V 



; in region II, (a) 



-i2k+u— ilk. 



-i2k+u—i2k—v. 



tanv ; in region III, {b) ( 9 ) 



in region IV, (c) 



where we have used two new separation constants k±, which are related to the previous ones 
k x and k y by the relation 4fc + /c_ = k 2 x + k 2 y . For convenience, we define also the following 
set of dimensionless constants: 



k-> 



L 1 L 2 k±, k 1 = JL 1 L 2 k 



LiL 2 h 



Even though S = {(u = 0, v < 0} U {u < 0, v 
scalar product is given by (see |L9f for details), 



(10) 

0} is a null hypersurface, a well defined 



— i / dxdy 



du + 



v=0 



L 9 v (j 



dv 



u=0 



'ir. 



Notice that the initial modes @ are well normalized on the boundary S between the flat 
region and the plane wave regions, and this means, from general grounds, that they remain 
well normalized on the boundary Si = {u = 0, < v < vr/2} U {0 < u < tt/2, v = 0} 
between the plane waves and the interaction region. Thus, the Cauchy problem for the 
interaction region is now well posed. 

However, although it has been rather easy to find the solution of the field equation in 
regions II and III which smoothly matches with the Minkowski vacuum, it turns out to be a 
difficult problem for the interaction region. Observe that the Cauchy data for the interaction 
region is imposed by the field modes (|9]a,b) on the lines S r = {u = 0, < v < 7r/2} U {0 < 
u < 7r/2, v = 0}, which are characteristic lines for the differential equation (||]). Thus, the 
only independent Cauchy data are the values of the function f(u,v) on them. Furthermore, 
we are only interested in finding the field solution on the causal past of the collision center, 
which is determined by the simple condition u = v. Then, the only relevant Cauchy data lie 

|27H for details). 



on the segments, Si = {u = 0, < v < ir/4} U {0 < u < 7r/4, 



0} 



sec 



In order to solve this partial problem, we start with the following change of coordinates, 



t = u + v, z = v — u, 



in equation (H) and we obtain, 



f, tt -f, zz + n(t,z)f = o, 

where the term Q(t, z), using (|I~0"D, is given by 

= ^ + (sin 2 t)/4 + k\ - (sin 2 z)/A 



cos 2 1 



COS 2 z 



'12) 



(13) 



(14) 



From now on, we will denote the causal past region of the collision center by S = {0 < t < 
7r/2, — 7r/4 < z < vr/4}. Observe that the behaviour of the variables t and z in equation 
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(|i~3|) in region S is very different. Since t runs from to ir/2 and z runs from —ir/A to 7r/4 in 
this region, the term flH| ) blows up as coordinate i goes to 7r/2, but is perfectly smooth over 
the entire range of coordinate z. This fact suggests that in the whole region S the physical 
results that we may expect are directly related to the coordinate t and we may not expect any 
physically remarkable change if we take z = in equation (|i3|). However, if we want to be 
consistent with such an approximation, we must also modify the boundary conditions that lie 
on the line segments Ei. Since on the boundary Ei we have that t = ±z and coordinate t runs 
from — 7r/4 to tt/4, we must also take t = z = 0. This means that the boundary conditions 
on Ei, given in (B|a,b), reduce in such an approximation to the flat boundary conditions 
given in (|9]c). Therefore we change the mode propagation problem for the colliding wave 
space-time into a rather simpler Schrodinger-type problem, which is clear from Fig. 2, and 
which requires only that we find a solution to equation ( |i"3"D with initial conditions given 
by the Minkowski flat modes below the hypersurface {t = 0, — 7r/4 < z < vr/4}. Recall, 
however, that none of the discussion above is applicable when a solution of equation ([13]) 
in the neighborhood of the folding singularities V and V is required. This is not only 
because in that case both coordinates t and z take values near 7r/2 and thus the potential 
term (|14]) is unbounded as z — > ir/2, but also because the boundary conditions (||a,b) are 
also unbounded as the folding singularities at t — ±z = it/ 2 are approached. In that case 
the mode propagation problem is much more complicated and a more detailed discussion is 



required (see |H5|, gg, ^ for details). 

In fact, rather than relying on the discussed approximations in the exact field equation 
(13), it will be necessary to rewrite a new field equation using an adequate approximation 
to the space-time geometry throughout the causal past region of the collision center. This 
is essentially because the process of renormalization involves the subtraction of the infinite 
divergences that arise from the formal definition of the stress-energy tensor, and these diver- 
gences can be expressed as entirely geometric terms, which are independent of any possible 
approximations in the field equation. 

Approximating the field equation (|T3"|) in the causal past of the collision center by taking 
z = is essentially equivalent to changing the line element, in the causal past of the collision 
center, by a related line element obtained from (f|) by setting z = 0, i.e., 

dsj = L1L2 (dt 2 — dz 2 ) — cos 2 tdx 2 — dy 2 . (15) 



We will suppose that the line element fll5|) applies all over the causal past of the collision 
center, not only in the interaction region but also through the plane wave regions II and III 
in the sense of Fig. 2. The plane wave collision starts at t = but to avoid smoothness 
problems derived from such an approximation, we will suppose that flTSp applies exactly on 
a range e < t < ir/2, for a certain e > 0. In the range < t < e, as described below, 
we will interpolate a line element which smoothly matches with the flat space at t = 0. 
Nevertheless, the details of this matching will not affect the main physical features. 
The exact field equation for this approximate space-time is, 

+ = 0, (16) 

where it is necessary to consider a coupling curvature term in the field equation because, 
although the exact space-time is a vacuum solution, we have a bounded nonzero value for R 
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in the approximated space-time. In order to solve this new field equation, we start rewriting 
the line element (|T5|) in the following general way, 

ds> = [hhh) df* - (M) dz 2 - dx 2 - (M) dy\ (17) 

where the /j are functions of coordinate t alone, which for values of < e < t < tt/2, 
can be straightforwardly determined by direct comparison with ( pT5[ ) as fi(t) = \[L\L 2 



f 2 (t) = \ / LiL 2 cost, fs(t) = cost. For values t < we take /i(t) = f 2 it) = 
/ 3 (t) = 1, which correspond to their values in flat space. Finally, in the interval < t < e, 
we smoothly interpolate each /j(t) (i = 1,2,3) between these values. Also, in order to 
prevent singularities in the field equation, we conveniently reparametrize coordinate t, by 
t*(t), as follows, 

dt* 1 



dt f 3 (ty K ' 

Now, we use the following ansatz for the field solutions, 

<p k = h(t*)e ikxX+ikvV+ikzZ , (19) 

where the plane wave factor in coordinates x, y is related to the translational symmetry of 
the space-time along the transversal directions x, y, and the plane wave factor in coordinate 
z is just a consequence of our approximation. Then equation fll6f) directly leads to the 
following Schrodinger-like differential equation for the function h(t*), 

Vt* + ^(t) h = 0, V(t) = co 2 (t) = f Q {t) + fl{t) k 2 x + f^it) k 2 y + flit) kl (20) 
where the function f Q (t) stands for, 

flit) = [Mt) Mt) f 3 it)} £R. (21) 

Such differential equation can be WKB solved, essentially because the short wavelength con- 
dition holds, i.e. uj~~ x d/dt* lncu 1. Observe that this condition reduces to (dt/dt*) dV/dt -C 
2u 3 , which becomes particularly accurate when the Killing-Cauchy horizon is approached 
since in that case dt/dt* = fsit) — > 0. Therefore, the mode solutions 0& which reduce to the 
flat mode solutions in the region prior to the arrival of the waves, are 

"1/2 . 
( j )k _ e ik x x+ik y y+ik 3 z~i J W(()dC (22) 

^i27T) 3 2k_Wit) 

where we denote u> 2 = k\ + k\ + k\ with k\ = \JL\Li k x , k 2 = \fLiL 2 k y , k 3 = k z and 
where W(t) stands for an adiabatic series in powers of the time- dependent frequency u{t) 
of the modes and its derivatives. Up to adiabatic order four (i.e. up to terms involving four 
derivatives of u(t)) Wit) it is given by, 

. . Ao Bo Aa Ba Ca Da 

Wit) =u+^ + ^r + ^r + ^ + ^ + -^, 23 
UJ A CU b UJ b CO 7 tU 9 LU 11 

where, using the notation V = dV/dt*, 
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4 32' 128 ' 258 ' 2048 ' 

and A n , B n , ... denote the n adiabatic terms in W(t). Up to adiabatic order zero it is simply 
W(t) = ui{t). Observe the two following facts: 

(i) Near the horizon t = tt/2 we have W(t) ~ u(t). This is because the higher adiabatic 
corrections vanish at the horizon. 

(ii) In the flat region prior to the arrival of the waves we have W(t) = Cj = (k\ + k\ + fcf) . 
In that case, since fz = 1, we can use (113) to set t* = t, where without loss of generality 



we choose the value t* — at t — 0. Therefore, the mode solutions (p2|) in the flat region 
reduce to, 

ilV ik x x+ik y y+ik z z—iu)t 

Yk ~ =(j 



(2tt) 3 2A;„ 

which indeed are the flat mode solutions defined in (|^c), recalling that the new separation 
constant k z = k 3 is related to the original k± by the ordinary null momentum relations, i.e., 

Cj = k + + k_, k z = k + — k_. (25) 

It is important to understand that we are constructing a set of mode solutions as an 
adiabatic series in terms of derivatives of the frequency u(t) in the differential equation (f2~0|). 
This procedure is similar but not equivalent to the construction of an adiabatic vacuum state 
where the field modes are expanded as an adiabatic series in terms of the derivatives of the 
metric coefficients (see for example [^j for details). In fact, observe for instance that the 



term /o(t) in (|20"D involves two derivatives of the metric since it is directly related to the 
curvature scalar. Thus, it would be a second order term for an eventual adiabatic vacuum, 
but it is simply a zeroth order term in our adiabatic series in derivatives of ui(t). 



4 Hadamard function in the interaction region 

The key ingredient to calculate the vacuum expectation value of the stress-energy tensor is 
the Hadamard function G^(x,x'), which is defined as the vacuum expectation value of the 
anticommutator of the field, i.e., 

G«(x, x') = ({(j>(x), 0(x')}> = E KO) <( x ') + u k{x') u* k (x)} . (26) 

k 

This Hadamard function contains non-physical divergence terms which can be subtracted 
by the following point-splitting prescription, 

G%\x,x') = G^\x,x') - S{x,x'), (27) 

where S(x, x') is the midpoint expansion of a locally constructed quantity commonly referred 
as a Hadamard elementary solution (see for example [|33[j ) and given by 
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S(x, x') 



2 
a 



2A< 2 > 



(2) 



a 



a 



af ln(/i-V) 



(o) 
;«2 o- 



ln(> V) - ^4 0) a 



(28) 



where the coefficients A^ 2 K P , A^ 



,(°) 



fiupri 



are written in Appendix A. We use the 



standard definition for the geodetic biscalar a(x,x') = (l/2)s 2 (x,x'), being s(x,x') the 
proper distance between the points x and x' on a non-null geodesic connecting them. Also, 
o~p,(x, x') = (dfdxP-)a(x, x') is a geodesic tangent vector at the point x with modulus s(x, x'), 
being x the midpoint between x and x' on the geodesic. The parameter fi in the logarith- 
mic term of (p8|) is an arbitrary length parameter, which is related to the two-parameter 
ambiguity of the point-splitting regularization scheme [Q. Then we can compute (T^) by 
means of the following differential operation, 



(T^(x)) = kD^^i'), 



(29) 



where T*^ is a nonlocal differential operator, which in the conformal coupling case (£ = 1/6) 
is given by, 



V 



[IV 



(30) 



However, the above differential operation and its limit have no immediate covariant meaning 
because G^'(x, x') is not an ordinary function but a biscalar and the differential operator 
T>^ v is nonlocal; thus we need to deal with the nonlocal formalism of bitensors (see, for 



29j for a review on this subject). 



example [^, [39| or the Appendix B of reference 

The regularization procedure (p7|), however, fails to give a covariantly conserved stress- 
energy tensor essentially because the locally constructed Hadamard function (^) is not in 
general symmetric on the endpoints x and x' (i.e. it satisfies the field equation at the point 



x but fails to satisfy it at x') (see |34| for details). Thus, to ensure covariant conservation, 
we must introduce an additional prescription: 



(TM) = (T*{x)) 



,(o), 



647T S 



(31) 



Note that this last term is responsible for the trace anomaly in the conformal coupling 
case, because even though (T^ u (x)) has null trace when £ = 1/6, the trace of {T^ v {x)) 

is given by (Tff) = —a2 > \x)/(lQn 2 ). The regularization prescription just given in ( |3~T|) 

57| , a set of properties that any 



satisfies the well known four Wald's axioms [|33|, |35l |36 



physically reasonable expectation value of the stress-energy tensor of a quantum field should 
satisfy. There is still an ambiguity in this prescription since two independent conserved local 
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curvature terms, which are quadratic in the curvature, can be added to this stress-energy 
tensor. In particular, the /i-parameter ambiguity in (ESI) is a consequence of this (see |B3 



for details). Such a two-parameter ambiguity, however, cannot be resolved within the limits 
of the semiclassical theory, it may be resolved in a complete quantum theory of gravity p3 |. 



Note, however, that in some sense this ambiguity does not affect the knowledge of the matter 
distribution because a tensor of this kind belongs properly to the left hand side of Einstein 
equations, i.e. to the geometry rather than to the matter distribution. 

After this preliminary introduction on the point-splitting regularization technique, we 
may proceed to calculate the Hadamard function G^\x, x') in the interaction region for the 
initial vacuum state defined by the modes <j>k, fl22|) . The Hadamard function can be written 
as, 

G^\x,x') = ^2M^4>l(x') +c.c. (32) 
k 

Note that solutions 0^ contain the function h(t*), which cannot be calculated analytically 
but may be approximated up to any adiabatic order as described in (p3|)-([2"4|). Thus, we have 
the inherent ambiguity of where to cut the adiabatic series. In fact, this is an asymptotic 
expansion, which has a well stablished ultraviolet limit but it may have convergence problems 
in the low-energy limit. However, observe from fll8|) and ( pOf) that since dt/dt* — > and 
V(t) = uJ 2 (t) — > k\ towards the horizon, the adiabatic series Q2"3|) reduces to W ~ u near 
the horizon. This means that we could cut the adiabatic series ( f23|) at order zero if we 
were interested in a calculation near the horizon. However, this is only partially true. In 
fact, it would be true if we were only interested in the particle production problem [|4j] but 



it is not sufficient for the calculation of the vacuum expectation value of the stress-energy 
tensor. This is because G^ calculated with h(t*) at order zero does not reproduce the 
short-distance singular structure of a Hadamard elementary solution (|28|) in the coincidence 
limit x — » x' . The smallest adiabatic order for the function h(t*) which we need to recover 
the singular structure of G^ 1 ' is order four, basically because our adiabatic construction of 
the mode solutions is similar (but not equivalent) to an adiabatic vacuum state (see |32| for 
details). 



Although expanding the function h{t*) in (|32| ) up to adiabatic order four will give an 
accurate value for the stress-energy tensor near the horizon, it will also give a suitable 
approximate value for this tensor all over the causal past of the collision center (region S in 
Fig. 2). The reason is that even though the short- wavelength condition, i.e. uj~ 1 d/dt* In uj <C 
1, is particularly accurate near the horizon it also holds throughout region S. 

In the mode sum ( |3~2"D we use the shortened notation Y^,k = Jo°° dk^/k^ f^dkx I^oo^ky 
or equivalently J2k = (-^i-^) -1 J^oo dk\ j x OD dk2 J^dk^/u, where the change of variables 
( |25| ) and the usual notation fllCf ) have been used. Therefore we have, 

G M(x,x') = * r r r dhdhdh x 



2(27r) 3 LiL 2 J-ooJ-ooJ-oo ^W(t)W(t' 
e _i /t*' W (O d C+ik x (x-x')+ik y (y-y')+ik z (z-z') _j_ c c ^33^ 

We assume that the points x and x' are connected by a non-null geodesic in such a way 
that they are at the same proper distance e from a third midpoint x. We parametrize the 
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geodesic by its proper distance r and with abuse of notation we denote the end points by 
x and x', which should not be confused with the third component of (t, z, x, y). Then we 
expand the integrand function in powers of e and we finally integrate term by term to get 
an expression up to e 2 . The details of such a tedious calculation can be found in ref. [j30| . 
The result is, 



G {1 \x,x') 



A + ab + C^ a V + D am a~ a ^^a l + -L { - \ - 2A^ pp ^ (34) 



-2A< 4 > 



(J,VpT 



a 



- ai 0) L 



in 2 1 o ^ a 



(2) 

1 )J,U 



(o) 



where L is a logarithmic term defined as L = 27 + ln(er £R/2), being 7 Euler's constant and 
where all the involved coefficients A, b, C^q... are given in Appendix A. 

According to (0), the Hadamard function can be regularized using the elementary 
Hadamard solution (123) and finally the regularized expression for G^\x,x') up to order 



e 2 is, 



1 



57T 



'3/9 

ai 0) L 



Q! „-/3 „-<5 



4 0) a( 2 ) 



+«P , „)»v-±<ip>a 



(35) 



where L is a bounded logarithmic term given by, L = 27 + ln(/i 2 £R/2), fi being the arbitrary 
length parameter introduced in (^), and where the coefficient B is B = b + 3c4 /(327r 2 ), 
which is given also in Appendix A. From (^) we can directly read off the regularized 
mean square field in the "in" vacuum state as ((f) 2 ) = A/2 — af^L/ilQix 2 ). It is important 
to remark, however, that the term in (^) appears only as a consequence of our 

approximate procedure of calculating the Hadamard function, i.e. using an adiabatic order 
four expansion for the initial modes in powers of the mode frequency uj(t) and its derivatives. 
Had we used an exact expression for the initial modes (or an adiabatic vacuum state 
such a term would not appear. 



5 Expectation value of the stress-energy tensor 

To calculate the vacuum expectation value of the stress-energy tensor we have to apply the 
differential operator (|30D to (^). As we have already pointed out, this is not straightforward 
because we work with nonlocal quantities. Note first that the operator fl3"U|) acts on bitensors 
which depend on the end points x and x', but the expression fl35|) for Gg depends on the 
midpoint x. This means that we need to covariantly expand (^) in terms of the endpoints 
x and x' . Also, the presence of quartic cx^ 1 terms in fl3"5|) gives, after differentiation, path 
dependent terms which must be conveniently averaged. The details of such a calculation may 
be found for instance in [^, ^] . Then, in the orthonormal basis Q\ = g\[ 2 dt, 62 = g x J 2 dz, 
63 = g]J 2 dx, #4 = g\jy dy, using the trace anomaly prescription (|31|), we obtain the following 
expectation values (T^ u ) in the conformal coupling case and for values 0<e<t<7r/2of 
coordinate t 
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<^> = 2880(W7r2 diag(-l,-l,l,-l) + -(r r )^ + 

diag( Pl (t), -p 2 (t), Pl (t) + 2p 2 (t), -p 2 (t)), (36) 

where we have used for simplicity the notation z = sint. The trace anomaly term in that 
case is given by ( ]3"I| ) as, 

l /TT v «2 _ diag(l, -1, -1, -1) 



4 X T '^ v 64 vr 2 ^ 2880(LiL 2 ) 2 tt 2 ' 

and the functions pi(t) and p 2 (t) are given by, 

-247247 + 1456234 z 2 + 792789 z 4 
PlW ~ 69189120 (LiL 2 ) 2 vr 2 (1 - z) 2 (l + zf ' ^ ' 

502931 + 179686 z 2 + 923887 ^ 4 
~ 69189120 (L X L 2 ) 2 vr 2 (1 - z) 2 (l + ^) 2 ' ( ' 

Recall that pi(t) is a positive definite function in an interval 0<e<t<7r/2. Both functions 
are unbounded at the horizon (t = vr/2), and the expectation value of the stress-energy tensor 
near the horizon is approximately given by, 

(?»U/2 = dia s 0(*), -p(t), Mt), -p(t)) + \{K)g^ (39) 

where p(t) = A(LiL 2 )~ 2 cost -4 and A ~ 0.0029. The behavior of (T M „) entirely agrees with 
the previous result |30| . For values of t < 0, (T^) = 0, and to be consistent with the 



approximation we have used for the space-time geometry, we should require that the value 
of {T^J) (0), which is valid for0<e<t<7r/2, goes smoothly to zero as t — > 0. In fact, 
this can be achieved using an adequate matching of the line element (|T5|) with the flat line 
element through the interval < t < e. 

Observe that the logarithmic term in the stress-energy tensor (|36D appears as a conse- 
quence of a similar term in the Hadamard function (p5|). The argument of this logarithm 
depends on the curvature scalar and thus it will grow unbounded as the flat region is ap- 
proached. However, the coefficient that will appear in front of such a logarithm in the 
Hadamard function (|35|) , depends only on locally constructed curvature terms. Therefore, 
with an adequate matching of the space-time geometry, this coefficient will also smoothly 
vanish towards the flat space region, below t — 0, and it will not give a contribution to 
the stress-energy tensor. The details of such a matching, however, will not affect the main 
features of the stress-energy tensor (0), particularly when the Killing-Cauchy horizon is 
approached. 

We must recall, however, that although the value (|36| ) for (T^) satisfies asymptotically 
the conservation equation near the Killing-Cauchy horizon, it does not satisfy exactly the 
conservation equation throughout region S, essentially because it is obtained by means of an 
approximation in the field modes. Nevertheless, we could obtain a truly conserved (T^), in 
the context of the present approximation, by solving the conservation equation considering a 
(T^J) given by the following set of non-null components {(T*(t)), (T*(t)), (T y (t))} 
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with the conditions: i) (T^(t)) = (T*(t)), which is compatible with (|36|) and it is a physical 
consequence of the isotropy of the metric ( ]T5| ) along the y-z directions, ii) trace anomaly 
condition, i.e. (T*(t)) = (Tfi - (T{(t)) -2{T z z {t)) } hi) the ansatz (T*(f)) = pi(t), which 
is the approximate value of (Tf(t)) obtained in our calculation. Finally, the conservation 
equation gives straightforwardly values for (T 2 2 (t)) and (T£(t)) which are compatible with 
the values (T*(t)) = P2(t) and (T*(t)) = —pi(t) — 2p 2 (t) obtained in our approximation. In 
particular, they have the same behaviour near the Killing-Cauchy horizon. 



Inspection of (|3g) shows that not only is the weak energy condition satisfied [4(J , which 



means that the energy density is nonnegative for any observer, but also the strong energy 
condition is satisfied. 



6 Conclusions 

We have calculated the expectation value of the stress-energy tensor of a massless scalar 
field in a space-time representing the head on collision of two electromagnetic plane waves 
throughout the causal past of the collision center and in the field state which corresponds 
to the physical vacuum state before the collision takes place. We have performed the cal- 
culations in this particular region essentially because, following the directions of a previous 
work [3~0|, we could introduce a suitable approximation to the space-time metric (see Fig. 



2). This approximation not only has allowed us to dramatically simplify the calculations 
but also to keep unchanged the main physical features, in particular the behavior of the 
stress-energy tensor near the Killing-Cauchy horizon of the interaction region. In fact, such 
an approximation is also valid for more generic plane wave space-times, and this will be the 
subject of a forthcoming paper. 

The results we have obtained are entirely compatible with the previous result 0] and 
they may be briefly described as follows: before the collision of the waves (T^ u ) = 0, which 
correspond to the lower edge Fig. 2. Then, after the collision the value of (T^ u ) starts 
to increase until it grows unbounded towards the Killing-Cauchy horizon of the interaction 
region. The weak energy condition is satisfied, the rest energy density is positive and diverges 
as cos -4 1. Two of the principal pressures are negative and of the same order of magnitude of 
the energy density. The strong energy condition is also satisfied, (T^) is finite but (T^) (T MI/ ) 



diverges at the horizon and we may use ref [41 on the stability of Cauchy horizons to argue 



that the horizon will acquire by backreaction a curvature singularity. Thus, contrary to 



simple plane waves, which do not polarize the vacuum |15|, [16|, the nonlinear collision of these 
waves polarize the vacuum and the focusing effect that the waves exert seems to produce, in 
general, an unbounded positive energy density at the focusing points. Therefore, when the 
colliding waves produce a Killing-Cauchy horizon, that horizon may be, in general, unstable 
by vacuum polarization. 

In the more generic case when the wave collision produces a spacelike singularity it seems 
clear that the vacuum expectation value of the stress-energy tensor will also grow unbounded 
near the singularity. In fact, in a forthcoming paper we will extend the approximation 
introduced in the present work to a more generic plane wave spacetime with the objective 
of more generally proving that the negative pressures associated to the quantum fields could 
not prevent the formation of the singularity. 
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A Useful adiabatic expansions 



The coefficients for the Hadamard function in (^), using for simplicity the notation z = sint, 
are: 



A 



B 



C, 



c 



mi 



D 

D- z 
D 



muni 



^~m) 



D 



1 




60LiL 2 7T 2 ' 




525 - 426 z 2 


+ 485 z A 


60480 {L x L 2 ) 2 t: 2 { 


1 -z) 2 {l + z) 2 


3597 + 2002 z 2 


+ 3107z 4 


332640 LiL 2 tt 2 (1 


-z) 2 (l + z) 2 


561 + 2222 z 2 


+ 1524 z A 



166320 (LiL 2 ) 2 7T 2 (1 -z)(l + z) 
1 



C s 



L1L2 

4433 + 12974 z 2 + 3513 z A 
823680 tt 2 (1-^) 2 (1 + z) 2 

5005 + 31200 z 2 + 45867 z 4 
4324320 {L l L 2 ) 2 "K 2 

1 



{L X L 2 ) 2 zzzz 

1001 + 4030 z 2 + 2848 z A 
864864LiL 2 7r 2 (l-^)(l + ^) 

143 + 260 z 2 + 120 z A 



i-xyy 



20592^^2 vr 2 (1 

D 7. 7 T V 



z) 2 1 + z) 



The coefficients for the midpoint expansion of the locally constructed Hadamard function 
D are: 
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Figure captions 

Fig. 1 The colliding wave space-time consists of two approaching waves, regions II and III, 
in a flat background, region IV, and an interaction region, region I. The two waves move 
in the direction of two null coordinates u and v. The four space-time regions are separated 
by the two null wave fronts u — and v = 0. The boundary between regions I and II is 
{0 < u < 7r/2, v = 0}, the boundary between regions I and III is {u = 0, < v < n/2}, and 
the boundary of regions II and III with region IV is S = {w < 0, v = 0} U {u = 0, v < 0}. 
Region I meets region IV only at the surface u — v — 0. The Killing- Cauchy horizon in the 
region I corresponds to the hypersurface u + v = n/2 and plane wave regions II and III meet 
such a Killing- Cauchy horizon only at V = {u = n/2, v = 0} and V' — {u — 0, v — n/2} 
respectively. The hypersurfaces u = n/2 in region II and v = n/2 in region III are a 
type of topological singularities commonly referred as folding singularities and they must be 
identified with V and V respectively. 

Fig. 2 The subset of Cauchy data which affects the evolution of the quantum field along 
the center u — v of the plane wave collision lies on the segments Si = {0 < u < n/A, v = 
0}U{u = 0, < v < n/A}. Region S is the causal future of this Cauchy data (or equivalently, 
the causal past of the collision center). 

We change the mode propagation problem for the plane wave collision, in the causal past 
of the collision center, region S, by a much simpler Schrodinger-type problem which consists 
in eliminating the dependence of the field equation on coordinate z by taking z = 0, and 
substituting the Cauchy data which come from the single plane wave regions on segments 
Si by much simpler Minkowski Cauchy data. This procedure is essentially equivalent to 
modifying the space-time geometry in the causal past of the collision center by eliminating 
the dependence on coordinate z in the line element, setting z = 0, and smoothly matching 
this line element, through plane wave regions II and III, with the flat spacetime below the 
segment {t = 0, — n/A < z < n/A}. 
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